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In the present work we calculate the theoretical tunneling conductance curves of SIN 
junctions involving high-Tc superconductors, for different possible symmetries of the 
order parameter (s, d, s + id, s + d, anisotropic s and extended s). To do so, we solve the 
real-axis Eliashberg equations in the case of an half-filled infinite band. We show that 
some of the peculiar characteristics of HTSC tunneling curves (dip and hump at eV > A, 
broadening of the gap peak, zero bias and so on) can be explained in the framework of 
the Migdal-Eliashberg theory. The theoretical dl /dV curves calculated for the different 
symmetries at T = 4 K are then compared to various experimental tunneling data 
obtained in optimally-doped BSCCO, TBCO, HBCO, LSCO and YBCO single crystals. 
To best fit the experimental data, the scattering by non-magnetic impurities has to be 
taken into account, thus limiting the sensitivity of this procedure in determining the 
exact gap symmetry of these materials. Finally, the effect of the temperature on the 
theoretical tunneling conductance is also discussed and the curves obtained at T = 2 K 
are compared to those given by the analytical continuation of the imaginary-axis solution. 

1 Introduction 

Both experimental evidences and theoretical arguments suggest that in high-Tc su- 
perconductors (HTSC) the /c-space symmetry of the energy|-|gap is different from 
the isotropic s-wave typical of conventional superconductors 111. A d-wave gap sym- 
metry with nodes on the Fermi surface, and some mixed symmetries, have been 
alternatively proposed as possible candidates for explaining the sometimes puzzling 
features of HTSC. 

The experimental determination of the symmetry of the energy gap is therefore 
a key step for understanding the nature of the pairing between electrons. In prin- 
ciple, the most usual and straightforward approach to this problem should consists 
in studying the tunneling conductance of HTSC junctions. 

In this paper, we calculate the theoretical tunneling conductance curves for 
superconductor - insulator - normal metal (SIN) junctions involving an high-Tc 
superconductor, for various symmetries of the energy gap. We then compare them to 
the experimental curves obtained on several high-Tc materials. For the calculation, 
we make use of the Migdal-Eliashberg theory □. The reasonabless of this approach 
is proven, among others, by the fact that, for suitable choices of the parameters, the 
calculations well reproduce several typical features of the experimental tunneling 
curves. 

2 The model 

Because of the layered structure of copper-oxide superconductors, we can restrict 
ourselves to a single Cu-0 layer, thus reducing the dimensionality of the problem. 
As a result, the fc-space becomes a plane and the Fermi surface is a nearly circular 



line on this plane. The quasiparticle wavevectors k and k' arc then completely 
determined by their modulus fcp and their azimuthal angles (fi and </>'. 

In order to calculate the tunneling conductance curves we solve the real-axis 
Eliashberg equations (EE) for the order parameter A(w, 0) and the renormaliza- 
tion function Z{uj,(j)). The kernels of these coupled integral equations contain the 
retarded electron,-|bQSon interaction a'^{Q, (p, (j)')F{Q) and the Coulomb pseudopo- 
tential /x* (0, (/)')I3'QB. To allow for various symmetries of the gap, we suppose that 
these quantities can be expanded, at the lowest order, in the following way: 

a^{n, cj,')F{^) = alF{n)iji, (0) t/;;, + al,F{n)^,^ (0) V^n (0') (1) 

fi* {(p, 0') = A^rsV-is (0) V-is (0') + A'InV'an (0) V-an {^') ■ (2) 

Here, the subscripts mean isotropic and anisotropic respectively, and V'is 
and i}}an {4>) are basis functions defined by : Vis (</>) = 1; V'an {4>) — cos (2(/)) for 
the d-wave, ?Aan (</') = 8-\/2/35 cos'* (2(/)) for the anisotropic s-wave, and V'an (0) = 
— 2-y/2/3 cos^ (20) for the extended s-waveEJ. 

We are interested in solutions of the real-axis EE containing separate isotropic 
and anisotropic terms, such as: 

A(t^, 0) ^ Ai,(tj) + Aa„(a;)V'an (</>) (3) 

Z{UJ, <p) = Zis{u}) + Zan(t^)'(/'an (0) • 

Inserting these expressions for A{uj,(j)) and Z{uj,(j)) in the real- axis EE makes 
them split into four equations for Ajs, Aan, ^is and Zan- From now on, we will put 
/i*=Q.for simplicity and we will admit that Zan is identically zero, as usually hap- 
pensu. The three remaining equations will be reported elsewhereu. To further sim- 
plify the problem we will put a^^F{fl) = g ■ a? F(17) where g is a constants. Then, 
Ais - 2/+°° dnalF{n)/n and Aan = (1/^) C d<^^an W Io°" dnal,F{n)/n re- 
suit to be proportional: Aan — 9 ■ Ms- 

Once obtained the real-axis solutions Ais(cLi) and Ai^^i'-^), we calculate the quasi- 
particle density of states: 



2tt 



Niu;) = ^ I dcjy (4) 

^Juj^- [Ais(a)) + Aan(t^)V'an (0)]' 



whose convolution with the Fermi distribution function gives the normalized con- 
ductance the main features of which will be discussed in the following section. 



3 Theoretical results 



The function af^F(fl), which is needed for numerically solving the real-axis EE, is 
actually unknown. As a first approximation, we use for afgF{fl) the experimental 
a^F {fl)-Bi22i2 (determined in a previous paper by starting from tunneling data on 
Bi2212 break junctionsB) suitably scaled to give Tc = 97 K (which is chosen as a 
representative value for HTSC): 

alF{n) = %a2F(17)Bi22i2. (5) 
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Figure 1: (a,b) Theoretical tunneling density of states for various symmetries at T=i K and ^*=0. 
The values of Ajg and Aan corresponding to each curve are indicated, (c) Theoretical normalized 
tunneling conductance in the d-wave symmetry, for increasing values of the coupling constants. 
The dip at about 2ypoak and the hump at about 3Vpeak are well evident. 



The details of the shape of a1i.F{il) are actually not relevant for the solutions, 
which instead strongly depend on the quantity wiog = exp /q^°° duj °' ^J^'^^ Incj^ . 

Figures |l|(a) and |l|(b) show the theoretical conductances calculated at the tem- 
perature r = 4 K for all the symmetries analyzed here: s, d, s + ic?, s + d, anisotropic 
s and extended s. For each curve, a couple of Ajg and Aan values is indicated which 
gives the expected symmetry and the correct critical temperature. As we shall see 
in the next section, these curves well reproduce some of the characteristic features 
of HTSC, such as the conductance excess below the gap and the broadening of the 
peak. Another apparently anomalous feature of HTSC recently reported in litera- 
turelZI is a suppression of the conductivity at about 2V^eak (usually known as "dip" ) 
and a subsequent enhancement at about SV^oak ( "hump" ) . Both these effects can be 
easily obtained by ^tarting from the Eliashberg equations, using high values of the 
coupling Constantsa and reducing uiog in order to maintain the critical temperature 
close to the experimental value. Figure ^(c) shows the conductance curves in the 
d-wave symmetry obtained for increasing values of the coupling constants. Note 
that the appearance of the "hump" is strictly related to the presence of the dip, as 
required by the conservation of the states. For instance, in the framework of the 
Eliashberg theory, a similar effect on the tunneling conductance curves jean also be 
obtained by simply taking into account the finiteness of the bandwidth B. 
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Figure 2: Comparison of the normalized tunneling conductance curves obtained by analytical 
continuation of the imaginary-axis solutions of the EE (dashed curves) to those obtained by direct 
solution of the real-axis EE (solid curves). All the curves are calculated at T=2 K. The agreement 
is good for the d, s -|- id and anisotropic s symmetries. In the s case, the analytical continuation 
shifts the gap of about 3 meV toward higher energies. In the s + d case the peak of the imaginary- 
axis solution is lower and broader than that of the real-axis one. The extended s curves agree quite 
well except that for very low voltage. 



The curves we have discussed so far were calculated at T—A K, in view of a 
comparison to the experimental data which are usually obtained at the temperature 
of liquid helium. We also calculated the N{V) curves at lower (T=2 K) and higher 
(r=40, 80 K) temperatures. The curves at T~2 K are not much different from those 
reported in Figures |^(a) and 0(b), apart from the fact that some fine structures (such 
as, for example, the first peak in the anisotropic s curve) are better distinguishable. 
Increasing the temperature results, in general, in a smoothing and broadening of 
the conductance peak at the gap edge. Above a temperature of about Tc/S all the 
conductances, apart the s-wave one, are almost the same. 

Just before comparing the theoretical results to the experimental data obtained 
on tunnel junctions, let us discuss an important point concerning the technique here 
adopted for the solution of the EE. The most usual and simplest approach to the 
solution of these equations consists in facts in solving them in their imaginary-axis 
formulation, and then to continue the solutions A(iwn) and Z(i(jJn) to the real axis. 
Actually, as we show in Figure |^, the results of this procedure only approximately 
agree with those obtained by directly solving the real-axis EE. Note that, in addi- 
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Figure 3: Comparison of the theoretical normalized tunneling conductance curves for the s + id- 
wave symmetry to those experimentally obtained on BSSCO (a), LSCO (b), YBCO (c), TBCO 
(d), HBCO (e). In each case the experimental and the calculated critical temperatures (Tc and 
j-calc^ are indicated. Here F is proportional to the impurity concentration and Co = cot((5o), where 
So is the scattering phase shift. 



tion, the analytical continuation is meaningful only at very low temperatures (this 
is the reason why the curves in Figure || are calculated at T=2 K). 

4 Comparison to experiments 

Let us now compare the-Lheoretical conductance curves to those obtained on optimally- 
doped single crystals^ of BSSCO , LSCO, YBCO, TBCO and HBCO . Here we 
only use for the comparison the theoretical curves obtained in the d symmetry, 
which in fact better agree with experimental data. For each of these materials aie 
use for afgF{n) the phonon spectral density determined by neutron scattering Eil, 
suitably scaled to obtain the best fit. As anticipated in a previous section, a large 
zero-bias often appears in the experimental tunneling curves. To best fit them, we 
add in the EE a term taking into account the scattering from impurities u. We then 
solve the real-axis EE with varying A^^ and Xam and choose the couple of values 
which give the best-fit curve. 

The results are shown in Figure ^ (a)-(e). In the case of BSSCO (a) and 
YBCO (c) the theoretical curves well agree with the experimental data, and also 
the calculated critical temperature is close to the measured one (Tc). In the case of 



LSCO (b) the observed dip above the gap cannot be reproduced in the framework 
of Ehashberg theory. In fact, the only way for obtaining the great zero-bias value 
(about 0.8) which appears in this dataset, is to increase the impurity concentration 
to such a large amount that the peak is broadened and the dip cancels out. For 
TBCO (d) and HBCO (e), the height of the peak at the gap edge is quite different 
in the experimental and theoretical curves. Actually, this arises from the choice of 
the normalization for the experimental data we took from the literature. With this 
normalization, in fact, the experimental curves seem to violate the conservation of 
the states, which is instead necessarily obeyed by the theoretical curves. In these 
cases, the agreement has to be evaluated by comparing only the shape of the curves, 
as they were plotted on different scales. Finally, while for TBCO the experimental 
and the theoretical critical temperatures coincide, in the case of HBCO we have 
j-caic^j^^^ K, which is much greater than the experimental one (Tc=97 K) 



5 Conclusions 



We discussed the effect of different possible symmetries of the order parameter 
on the tunneling conductance curves of HTSC in the framework of the Eliashberg 
theory. We showed that this theory allows to explain many characteristic features 
of HTSC which differ from a classic BCS behaviour, and is therefore a useful tool 
to investigate the properties of these materials. We also tested the theoretical 
expectations by comparison with experimental data obtained on many different 
copper-oxide superconductors, and we found in many cases a reasonable agreement 
between the results of calculation and measurements. The most probable symmetry 
of the order parameter seems to be the c?-wave, even though the sensitivity of this 
method in distinguish between different possible symmetries of the pairing state is 
strongly limited by the quality of the crystals. 
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